Directed and elliptic flow in heavy ion collisions from i?beam = 90 Me V/nucleon to 

E c , m , = 200 GeV/nucleon 
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Recent data from the NA49 experiment on directed and elliptic flow for Pb+Pb reactions at 
CERN-SPS are compared to calculations with a hadron-string transport model, the Ultra-relativistic 
Quantum Molecular Dynamics (UrQMD) model. The rapidity and transverse momentum depen- 
dence of the directed and elliptic flow, i.e. t>i and «2, are investigated. The flow results are com- 
pared to data at three different centrality bins. Generally, a reasonable agreement between the 
data and the calculations is found. Furthermore, the energy excitation functions of Wi and V2 from 
-Ebeam = 90A MeV to -Bcm = 200A GeV are explored within the UrQMD framework and discussed 
in the context of the available data. It is found that, in the energy regime below -Ebcam < WA GeV, 
the inclusion of nuclear potentials is necessary to describe the data. Above 40A GeV beam energy, 
the UrQMD model starts to underestimate the elliptic flow. Around the same energy the slope of 
the rapidity spectra of the proton directed flow develops negative values. This effect is known as 
the third flow component ("antiflow") and cannot be reproduced by the transport model. These 
differences can possibly be explained by assuming a phase transition from hadron gas to quark gluon 
plasma at about 40A GeV. 
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I. INTRODUCTION 

To create extremely hot and dense matter with par- 
tons as its fundamental components - called the Quark- 
Gluon Plasma (QGP) - is a major goal of current and fu- 
ture high energy heavy-ion collisions experiments at SPS, 
RHIC and LHC [jj. However, due to the complex nature 
of the relativistic nucleus-nucleus reactions, the QGP, if 
it has been created, escapes direct detection. Therefore, 
in order to establish the existence and later on to in- 
vestigate the properties of the new state of matter, one 
must find observables which allow to deduce the proper- 
ties of the intermediate (QGP) state from the final state 
hadrons. 

The exploration of the transverse collective flow is the 
earliest predicted observable to probe heated and com- 
pressed nuclear matter Q. The transverse flow is inti- 
mately connected to the pressure gradients. Therefore, 
it is sensitive to the equation of state (EoS) and might 
be used to search for abnormal matter states and phase 
transitions dUl]. 

The intermediate energy regime available at CERN- 
SPS or at the future GSPFAIR facility is often referred 
to as the right place to look for a phase transition to the 
QGP. Lattice QCD (1QCD) calculations @, Q show that 
the critical temperature is around 170 MeV (for fij, = 
0) and the critical energy density is around 1 GeV/fm 3 . 
These values can already be reached at 20-30 AGeV beam 
energy. At finite baryo-chemical potential, the heated 
and compressed nuclear matter created at these energies 
crosses the phase transition line possibly even on the high 
fi side of the critical endpoint. Therefore, it is possible to 



talk about a phase transition of first order, here. During 
such a first order phase transition the softest point in 
the equation of state would be mostly pronounced. For 
example, the proton antiflow around midrapidity ( "third 
flow component" Q ) and the collapse of the elliptic flow 
observable have been declared as a signal for the phase 
transition (9l. Hoj. 

In this paper, recent results on proton and pion di- 
rected and elliptic flow from the NA49 experiment [llj 
are investigated and predictions for FAIR are presented. 
The data is compared to transport model calculations 
(UrQMD v2.2). The proton flow measures the behaviour 
of the nuclear matter during a heavy ion collision, while 
the flow of pions is a sign for the properties of newly 
produced particles. 

The paper is organized as follows. Section II includes 
an introduction of the UrQMD model. Section III and 
IV introduce the flow systematics and the different mea- 
surement methods. Then, in Section V, the directed flow 
results are shown. There are rapidity and transverse mo- 
mentum distributions for 40 and 160 AGeV beam energy. 
The centrality dependence is also studied. Predictions for 
rapidity and p t dependence of v\ for Eb eam = 20A GeV 
and 30^4 GeV are made. In section V (C) the energy de- 
pendence of the slope around midrapidity of the directed 
flow is investigated in the context of the available data. 
Afterwards, in Section VI, the same analysis for the ra- 
pidity and transverse momentum dependence of elliptic 
flow (1*2) is shown. Section VI (C) discusses the excita- 
tion function of elliptic flow over the whole energy range 
from SIS to RHIC. Section VII summarizes the paper. 
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II. THE URQMD MODEL 

For our investigation, the Ultra-relativistic Quantum 
Molecular Dynamics model (UrQMD v2.2) [H, El i s a P~ 
plied to heavy ion reactions from £"beam = 90A MeV 
to ^snn = 200 GeV. This microscopic transport ap- 
proach is based on the covariant propagation of con- 
stituent quarks and diquarks accompanied by mesonic 
and baryonic degrees of freedom. It simulates multiple 
interactions of in-going and newly produced particles, the 
excitation and fragmentation of colour strings and the 
formation and decay of hadronic resonances. Towards 
higher energies, the treatment of sub-hadronic degrees of 
freedom is of major importance. In the present model, 
these degrees of freedom enter via the introduction of a 
formation time for hadrons produced in the fragmenta- 
tion of strings [3, HE Ell • A phase transition to a quark- 
gluon state is not incorporated explicitly into the model 
dynamics. However, a detailed analysis of the model in 
equilibrium, yields an effective equation of state of Hage- 
dorn type [llGJ]. 

The UrQMD transport model is successful in describ- 
ing the yields and the pt spectra of different particles 
in pp and pA collisions [19j. It has also been ap- 
plied to study the flow at lower energies and at RHIC 
energies ( [U M M M M HI ) . 

III. FLOW SYSTEM ATICS 

The first coefficient of the Fourier expansion of the 
azimuthal distribution of the emitted particles (v±) de- 
scribes the directed in-plane flow. The directed flow mea- 
sures the total amount of transverse flow. It is most pro- 
nounced in semi-central interactions around target and 
projectile rapidities where the spectators are deflected 
away from the beam axis due to a bounce-off from the 
compressed and heated matter in the overlap region, vi 
is defined by 

V t = (COS(0 - $kp)) , (1) 

where <fi denotes the azimuthal angle of one outgoing par- 
ticles and Qrp is the azimuthal angle of the reaction 
plane. The angular brackets denote an average over all 
considered particles from all events. 

Three different interesting properties of the directed 
flow have been proposed. 

• The time scales probed by the directed flow are set 
by the crossing time of the Lorentz-contracted nu- 
clei. Thus, it serves as keyhole to the initial, prob- 
ably non-equilibrium, stage of the reaction [261 ] . 

• The softening of the equation of state in a first order 
phase transition leads to a decreasing directed flow 

m, M Si- 

• The space-momentum correlation of the emitted 
particles can be adressed experimentally via the v\ 
rapidity distributions of nucleons and pions. 



The second coefficient of the Fourier expansion of the 
azimuthal distribution of the emitted particles (v^) is 
called elliptic flow H, M,WH H M H HE H GH S [$1 . 
This type of flow is strongest around central rapidi- 
ties in semi-peripheral collisions. It is driven by the 
anisotropy of the pressure gradients, due to the geomet- 
ric anisotropy of the initial overlapping region. There- 
fore, it is a valuable tool to gain insight into the expand- 
ing stage of the fire ball. v 2 is defined by 

v 2 = (cos[2(0-$ RP )]) . (2) 

There are two competing effects which lead to contri- 
butions with different signs to the integrated vi value. At 
low energies or early times there is the so called "squeeze- 
out" -effect. The spectator matter blocks the emission in 
the impact parameter direction and therefore the flowing 
matter is "squeezed" -out perpendicularly to the reaction 
plane. This leads to negative elliptic flow values. The 
second effect is the so called in-plane flow. This type of 
flow becomes important at higher energies and/or later 
times. At higher bombarding energies (£a a b > 10^4 GeV) 
the spectators leave the interaction zone quickly. The 
remaining hot and dense matter expands almost freely, 
where the surface is such that in-plane emission is pre- 
ferred. Therefore the elliptic flow receives a positive con- 
tribution. 

Let us now explore the time evolution of the pressure 
gradients in connection with the elliptic flow develop- 
ment. The transverse pressure gradients have been cal- 
culated for the first 10 fm at E iab = AOA GeV (see Fig. [T]) 
and the highest SPS energy (see Fig. [2]). In both cases 
one observes large pressure gradients in the very early 
stage of the collision. For the lower energy the maximum 
is reached around t = 3 fm and for the higher energy it 
is shifted to even earlier times. The difference between 
the pressure gradients in x- and y- direction is respon- 
sible for the vi development. As it can be seen in Figs. 
QJbottom) and [2fbottom) the temporal evolution of el- 
liptic flow for pions starts exactly after this maximum. 
The elliptic flow increases during <~ 6 fm until it reaches 
almost its final value. After t — 10 fm it decreases a little 
because of resonance decays. So, elliptic flow builds up 
in the early stage of the collision due to the difference of 
pressure gradients as it is expected. 

IV. FLOW MEASUREMENT METHODS 

In the UrQMD model calculation of flow, the exact az- 
imuthal angle of the reaction plane $rp is known by defi- 
nition. However, an unambiguous experimental measure- 
ment of the azimuthal anisotropic flow is not a trivial task 
due to the unknown orientation of the reaction plane. Of- 
ten, experiments use the so called reaction plane method 
[4fj| to extract the magnitude of flow. In this method, 
the reaction plane is fixed according to the flow vector 
of the event, then the estimated v 2 with respect to the 
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FIG. 1: (Color online) UrQMD calculation for the time evolu- 
tion of the pressure gradients and elliptic flow for Pb+Pb 
interactions at Ei^b = 40 A GeV. Top: dP/dx (full line), 
dP/dy (dotted line) and the difference between these two AP 
(dashed line) are depicted. Bottom: Elliptic flow of pions 
(squares) versus time at midrapidity for mid-central collisions 
(b=7fm). 



FIG. 2: (Color online)UrQMD calculation for the time evo- 
lution of the pressure gradients and elliptic flow for Pb+Pb 
interactions at _Ei a b = 160 A GeV. Top: dP/dx (full line), 
dP/dy (dotted line) and the difference between these two AP 
(dashed line) are depicted. Bottom: Elliptic flow of pions 
(squares) versus time at midrapidity for mid-central collisions 
(b=7fm). 



chosen reaction plane is corrected for the event plane res- 
olution, which accounts for the error in the deduction of 
the reaction plane. 

However, these two-particle correlations based method 
might suffer from effects which are not related to the re- 
action plane, these additional contributions are usually 
called non-flow effects [4l| , such as the overall transverse 



momentum conservation, small angle azimuthal correla- 
tions due to final state interactions, resonance decays, 
jet production [i^] and quantum correlations due to the 
HBT effec t 1431 ] . Recently, the cumulant method was pro- 
posed 0, l45l] to diminish the non-flow effects. The idea 
of the cumulant method is to extract flow with many- 
particle cumulants, which are the many-particle corre- 
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lations with subtraction of the contributions from the 
correlations due to the lower-order multiplets. It is be- 
lieved that the pure many-particle non-flow correlations 
have much less contributions to the measured flow in 
the many-particle cumulant method. In other words, the 
many-particle cumulant method should be much less sen- 
sitive to non-flow effects 0, EH • 

This has been confirmed by the test of the cumulant 
method with the UrQMD model [U[25j]. However, the 
test also shows that, at least for the V2 measurement, the 
cumulant method is not completely free from the effect 
of event-by-event V2 fluctuations [43]. Especially when 
the genuine v-i signal is weak (for example, for the most 
central events and very peripheral events), the effect of 
V2 fluctuations could be so strong that the results from 
the cumulant method becomes unreliable. 



V. DIRECTED FLOW RESULTS 
A. Rapidity dependence 

Fig. [3] shows the rapidity dependence of the directed 
flow of protons for central (b < 3.4 fm), mid-central (6 = 
5 — 9 fm) and peripheral (b > 9 fm) Pb+Pb reactions 
at £q a b = 40 A GeV. The symbols denote data by the 
NA49 collaboration analyzed with different methods [ll[ . 
UrQMD calculations are depicted with black lines. The 
directed flow is most pronounced at high rapidity values 
where the bounced off spectator matter sits. 

The strong centrality dependence of v% can be seen, 
as the shape of the curves changes drastically from cen- 
tral to peripheral collisions. Overall, the model calcu- 
lations are in line with the data from the event plane 
method (standard method, full circle). It is interesting 
to note that one kind of non-flow effects, i.e. momentum 
conservation, has been subtracted in this reaction plane 
method already [48| . Unfortunately, the 2-particle cumu- 
lant measurements seem to be affected by the non-flow 
effects which are in this analysis dominated by momen- 
tum conservation. This effect is most pronounced for 
peripheral collisions. Here, the experimental data points 
do not approach zero at midrapidity For the more reli- 
able 3-particle cumulant method [49(, the experimental 
results agree well with those of the reaction plane method 
within the statistical error (please see also the following 
v\ data). This indicates that the momentum conserva- 
tion corrected reaction plane method gives also reliable 
V\ data. 

Fig. [4]shows directed flow of pions at E\ a b — 4QA GcV. 
The different sign with respect to the proton flow can be 
explained by shadowing. The pions are newly produced 
mesons and therefore they are composed by a quark and 
especially an anti-quark. Therefore, the cross-section of 
the pions with the nuclear matter is so large that they 
cannot escape in the direction where the rest of the col- 
liding nuclei/spectator matter resides. 

For peripheral collisions the measured directed flow 
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FIG. 3: (Color online) Directed flow of protons in Pb+Pb col- 
lisions at £i ab = 40 A GeV withp t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 



reaches the same amount as for protons (about 10% at 
y w 2), but for central collisions it is much less (only 
about 2.5% at y « 2). The UrQMD calculations overes- 
timate the pion directed flow at large rapidities by about 
a factor of two for all centralities. This overestimation 
might be explainable if one assumes that the NA49 col- 
laboration is not able to measure all produced pions. Es- 
pecially pions produced from nucleons which fly straight 
ahead through the collision producing only one or two pi- 
ons appear only in the veto calorimeter. However, in the 
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FIG. 4: (Color online)Directed flow of pions in Pb+Pb colli- 
sions at £i ab = 40 A GeV with p t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 
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FIG. 5: (Color online) Directed flow of protons in Pb+Pb col- 
lisions at £ lab = 160^1 GeV with p t < 2 GeV/c. UrQMD cal- 
culations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 
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model calculations every produced pion within the given 
rapidity and transverse momentum bin is taken into ac- 
count. 

In Fig. directed flow of protons at — 160A GeV 
is shown. In this case the model calculations slightly 
overestimate the flow at higher rapidities in contrast to 
the proton flow at -Ei ab = 40A GeV. Surprisingly, the 
data stays almost constant at zero. Even in peripheral 
collisions there are only about v\ — 2.5% at y — 1.5 for 
the reaction plane method analysis. The flat shape of 



the curve with respect to the lower energy results can be 
explained by the spectators sit at higher rapidity values 
in this plot. This spectator matter is responsible for the 
directed flow near beam/target rapidities. Since there 
are no data points of the reaction plane method above 
y = 1.5 this increase is not seen. 

For the pion directed flow at the higher SPS energy 
(Fig. [5]), one observes a reasonable agreement between 
the model calculations and the data from all three mea- 
surement methods, i.e. reaction plane method, i>i{2} 
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FIG. 6: (Color online)Directed flow of pions in Pb+Pb colli- 
sions at _Bi ab = 160A GeV with p t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



and i>i{3}. Above y = 1, the UrQMD results increase 
slightly stronger than the data. The centrality depen- 
dence is reproduced correctly. Both in the model and 
in the experimental data there is an increase in flow for 
more peripheral collisions. 

Fig. [7] and Fig. [5] are predictions for the rapidity 
dependence of proton and pion directed flow at E\ a \, — 
2QA GeV and E uh = 30A GeV. The shown results are 
calculated using the UrQMD model with a transverse 
momentum cut of pt < 2 GeV/c. In the present model, 
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FIG. 7: Directed flow of protons in Pb+Pb collisions at i?i a b = 
2QA GeV and £ lab = 30A GeV with p t < 2 GeV/c. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b — 5 — 9 fm for mid- 
central and b — 9 — 15 fm for peripheral collisions have been 
used. 



the calculated directed flow results at these energies un- 
der investigation at the new FAIR facility at GSI look 
rather similar to that at E\ a ^ = 4QA GeV. There is also 
an inverse centrality dependence for protons and pions 
in the amount of the directed flow. It is very interesting 
to see if the experimental proton flow data at this en- 
ergy will show the negative slope around midrapidity as 
predicted. 



B. Transverse momentum dependence 

Fig. [5] shows directed flow of protons at Ei a b = 
A0A GeV as a function of transverse momentum. For 
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FIG. 8: Directed flow of pions in Pb+Pb collisions at _E"i a b = 
20 A GeV and £ lab = 30 A GeV with p t < 2 GeV/c. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b — 5 — 9 fm for mid- 
central and b — 9 — 15 fm for peripheral collisions have been 
used. 



the calculation the same rapidity cut (0.1 < y < 1.1) as 
in the data obtained with the cumulant method has been 
used. For the standard reaction plane analysis the differ- 
ent cut of 0.1 < y < 1.8 has been applied. The different 
cuts have been used by the experimental collaboration 
to improve the statistics and reduce the time amount for 
the analysis. Since the higher order cumulant method 
measurements are the most reliable way to reduce sys- 
tematic errors 0, [HI, the comparisons are done with 
the cuts for the cumulant method. The calculations are 
in line with the fi{3} data for mid-central and periph- 
eral collisions. Unfortunately the statistical errors are 
the biggest for this method and there are only a few data 
points because it is the newest and most elaborate anal- 



0.2 
0.15 

0.1 
0.05 

0.0 
-0.05 



central 



i | i | i | i | i | i r 
Protons, 40 AGeV, 0.1 < y < 1.1 



UrQMD-2.2 

• standard method. 

- v,{2} 




0.6 0.8 1.0 1.2 1.4 
p t [GeV/c] 



2.0 



FIG. 9: (Color online) Directed flow of protons in Pb+Pb col- 
lisions at _Bi ab = 40A GeV with 0.1 < y < 1.1. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 



ysis. 

The transverse momentum dependence of the pion di- 
rected flow at i?iab = 40 A GeV (Fig. [TU]) shows many 
uncertainties. The data differs very much depending on 
the analysis method. Furthermore, the cumulant method 
of order 3 has large statistical error bars. The UrQMD 
calculations are negative as the sign convention suggests. 
It is remarkable that the directed flow does not show any 
clear transverse momentum dependence. 

Next, we turn to the results for 16QA GeV colli- 
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FIG. 10: (Color online) Directed flow of pions in Pb+Pb colli- 
sions at S lab = 40 A GeV with 0.1 < y < 1.1. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 
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FIG. 11: (Color online)Directed flow of protons in Pb+Pb 
collisions at £ lab = 160A GeV with 0.1 < y < 1.6. UrQMD 
calculations are depicted with black lines. The symbols are 
NA49 data from different analysis methods. The standard 
method (circles), cumulant method of order 2 (squares) and 
cumulant method of order 3 (triangles) are depicted. The 
12.5% most central collisions are labeled as central, the cen- 
trality 12.5% -33.5% as mid-central and 33.5% -100% as pe- 
ripheral. For the model calculations the corresponding impact 
parameters of b < 3.4 fm for central, b — 5 — 9 fm for mid- 
central and b — 9 — 15 fm for peripheral collisions have been 
used. 



sions. For central collisions the calculated transverse 
momentum dependence of the directed flow of protons 
at Eiah — 160A GeV (see Fig. [11]) starts at zero for 
Pt = and increases steadily until 2% at p t — 2 GeV/c. 
The measured data fluctuates between +2% and —2%. 
Going to mid-central collisions the calculations are quite 
in line with the third order cumulant method and the 
standard method data. In peripheral collisions there are 
big systematic uncertainties because especially around 
Pt = 1.8 GeV/c the data points differ between —4% and 



+5.5%. The second order cumulant results are higher, 
but they are not corrected due to momentum conserva- 
tion. The negative directed flow of protons is consistent 
with the rapidity distribution. Because of the negative 
slope around midrapidity at the higher SPS energy the 
protons fly in the "wrong" direction. 

The calculated directed flow for pions at Bi a b = 
160A GeV, as depicted in Fig. [T^l stays constant as a 
function of the transverse momentum at (—1)% for cen- 
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FIG. 12: (Color online) Directed flow of pions in Pb+Pb colli- 
sions at £ lab = 160A GeV with 0.1 < y < 1.6. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, 6 = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 
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FIG. 13: Directed flow of protons in Pb+Pb collisions at 
E la . h = 20 A GeV and £i ab = 30 A GeV with 0.1 < y < 1. 
UrQMD calculations for 20 AGeV are depicted with solid lines 
while the results for 30 AGeV are depicted by dashed lines. 
Impact parameters of b < 3.4 fm for central, b = 5 — 9 fm for 
mid-central and b = 9 — 15 fm for peripheral collisions have 
been used. 
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tral, at (—2)% for mid-central and at (—3)% for periph- 
eral collisions. In mid-central and peripheral collisions 
this trend is in line with the third order cumulant mea- 
surement. Especially the data from second order cumu- 
lant increase up to 4% — 6%. This is probably due to 
momentum conservation which is not taken into account 
to correct the data. The negative directed flow of pions 
is consistent with the rapidity dependence in the chosen 
bin. 

Finally, we show the prediction for the transverse mo- 



mentum dependence of directed flow at E\ a b = 20 A GeV 
and E U b = 30 A GeV in Figs. [13] and [H In the present 
calculations, the shape and magnitude of the flow are 
similar to the results at 4QA GeV. The directed flow of 
pions at 2QA GeV and 30yl GeV looks also rather similar 
to the calculations for the pion directed flow at AQA GeV. 
But there is a difference for peripheral collisions. The 
v i(Pt) value is only about (—2.5)% compared to (—5)% 
at E Uh = 40A GeV. 
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FIG. 14: Directed flow of pions in Pb+Pb collisions at £"i a b = 
20 A GeV and £ lab = 30 A GeV with 0.1 < y < 1. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b = 5 — 9 fm for mid- 
central and b = 9 — 15 fm for peripheral collisions have been 
used. 
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FIG. 15: (Color online)Slope of «i(t/) of protons around 
midrapidity extracted from normalized (y = y/yb) rapidity 
distributions. The data are taken from E895 (squares) [50] 
and NA49 (triangles) [ll[ . UrQMD calculations with included 
mean field (HMw)are depicted with full circles. Open circles 
depict UrQMD calculation in the cascade mode. 



the data. Here is shown the calculation with included 
mean field from a hard equation of state with momen- 
tum dependence and medium-modified nucleon-nucleon 
cross sections (HMw) [5l|, [HJ . At higher energies the cal- 
culation has been performed in the cacsade mode without 
the additional inclusion of nuclear potentials. 



C. Excitation function 

To characterize the amount and the direction of the 
directed flow of protons over the energy range from 
2 — 160 AGeV one can extract the slope around midrapid- 
ity from the normalized rapidity distributions usually re- 
ferred to as the "F" parameter [5(| ■ "Normalized" means 
in this case y/yb where yb is the beam rapidity. This nor- 
malization accounts for the trivial energy dependence of 
the slope. The values for the slope in Fig. [15] have been 
extracted via a polynomial fit of the form a ■ x + b ■ x 3 
with x = y/yb- At low energies, one observes that the 
inclusion of a nuclear potential is needed to reproduce 



At SPS energies the data develops even negative values 
for the slope around midrapidity 0, U3] • This behaviour 
cannot be reproduced within the transport model calcu- 
lation. However, ideal hydro calculations have predicted 
the appearance of a so-called "third flow component" [|| 
or "antiflow" [46| at finite impact parameters. In these 
analysis' it was pointed out that this "antiflow" devel- 
ops if the matter undergoes a first order phase transition 
to the QGP. In contrast, a hadronic EoS without QGP 
phase transition did not yield such an exotic "antiflow" 
(negative slope) wiggle in the proton flow V\(y) at low 
energies. 



11 



> 



central 



> 



mid-central 




UrQMD-2.2 
standard method - 

v 2 {2) 



> 



0.1 
0.08 
0.06 
0.04 
0.02 

0.0 

-0.02 

-0.04 

-0.06 

-0.08 

-0.1 
0.1 

0.08 

0.06 

0.04 

0.02 

0.0 

-0.02 

-0.04 

-0.06 

-0.08 

-0.1 
0.1 

0.08 

0.06 

0.04 

0.02 

0.0 
-0.02 
-0.04 
-0.06 
-0.08 

-0.1 



FIG. 16: (Color online)Elliptic flow of protons in Pb+Pb col- 
lisions at £ lab = 40A GeV withp t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, b = 5 — 9 frn for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 



VI. ELLIPTIC FLOW 
A. Rapidity dependence 

Elliptic flow develops because of the almond shape of 
the overlapping region in a heavy ion collision. The 
breakdown of proton elliptic flow at E\ a b = 40A GeV 
has been stressed as a signal for the observation of a first 
order phase transition [10]. As can be seen from the 




NA49 data in Fig. [T6l the elliptic flow parameter van- 
ishes at midrapidity only for the standard reaction plane 
method data, which could be affected by the non-flow 
effects. The cumulant measurements show a completely 
different shape. The second and fourth order vi measure- 
ments for peripheral and mid-central collisions increase 
to about 6 — 8% around midrapidity not consistent with 
zero. Looking now at the UrQMD calculations one can- 
not observe an overestimation, but finds that the results 
are compatible with the data. Therefore, before one can 
draw concise conclusions the systematic and statistical 
uncertainties of the data must be resolved. 

For the elliptic flow of pions (Fig. [TT]) in central colli- 
sions there are again systematic uncertainties in the ex- 
perimental data. Again all three methods have differ- 
ent shapes and the calculated elliptic flow value is con- 
sistent with the measurement with errors. The elliptic 
flow increases to about 1% around midrapidity within the 
UrQMD model. For mid-central collisions the picture be- 
comes clearer. The data and the calculations are in line. 
Here, the elliptic flow around midrapidity is about 2—3%. 
Going to peripheral collisions the experimental results of 
the standard method show an unsteady behaviour but 
are compatible with the UrQMD calculations. 

The elliptic flow is expected to be larger in more pe- 
ripheral collisions because the anisotropy in coordinate 
space which is the source of this flow component is larger. 
This dependence can be seen very well in the proton flow 
at i?i a b = 160A GeV (Fig [T5|) . In central collisions there 
is almost no elliptic flow at all. Going to mid-central 
collisions the flow is between 2 — 6% depending on ra- 
pidity and in peripheral collisions even between 2 — 8%. 
The model calculations do not show a big difference be- 
tween mid-central and peripheral collisions. Both values 
are about 2% around midrapidity. The large elliptic flow 
data in peripheral collisions could be due to the fact of 
non-flow effects and/or large event-by-event v-i fluctua- 
tions there, which affect the vi measurement heavily. 

Fig. 033 shows the elliptic flow of pions at = 
160A GeV. These newly produced particles show almost 
no elliptic flow in central collisions, i.e. they are emitted 
isotropically out of the collision region. The experimen- 
tal result is between 1 — 2%. In mid-central collisions 
there is a very good agreement between the model and 
data again. For peripheral collisions the same under- 
estimation of the calculations as for the proton flow in 
comparison to the data can be seen. 

Fig. [20] and Fig. [5T] are predictions for the rapidity 
dependence of proton and pion elliptic flow at £q a b = 
20A GeV and E lah = 30A GeV. The shown results are 
calculated using the UrQMD model with a transverse 
momentum cut of pt < 2 GeV/c. The elliptic flow results 
at this energy of the new accelerator FAIR at the GSI 
look rather similar to that at E\ a b — 40/1 GeV. 
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FIG. 17: (Color online)Elliptic flow of pions in Pb+Pb colli- 
sions at £i ab = 40 A GeV with p t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b — 9 — 15 fm for peripheral collisions have been used. 
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FIG. 18: (Color online)Elliptic flow of protons in Pb+Pb colli- 
sions at £ lab = l&QA GeV with p t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



B. Transverse momentum dependence 

The elliptic flow of protons at E iah = 160A GeV (Fig. 
|2"4"| in central collisions looks similar to that at A0A GeV. 
There is a smooth increase from zero at pt = to about 
2% at pt — 2 GeV/c. The underestimation of the flow 
by the transport model calculations at higher energies 
is visible in the results for mid-central and peripheral 
collisions at high p t . At low pt the calculations are in 
line with the data. 



For the elliptic flow of protons as a function of trans- 
verse momentum there is a different rapidity cut used 
in the data than in the calculations. For proton and 
pion flow at E\ ah = A0A GeV it is -0.1 < y < 1.1 
and for E lab = 160A GeV it is 0.1 < y < 1.6 for the 
cumulant order measurements. For the standard reac- 
tion plane method data, larger rapidity bins have been 
used for the integration. To improve the statistics of 
the UrQMD results we use symmetric rapidity cuts as 
\y\ < 1.1 for 40A GeV and \y\ < 1.6 for 160A GeV. 
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FIG. 19: (Color online)Elliptic flow of pions in Pb+Pb colli- 
sions at _Bi ab = 160A GeV with p t < 2 GeV/c. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



Since the elliptic flow is symmetric in rapidity this can 
be done without problems. The elliptic flow of protons 
at £i a b = 40A GeV (Fig. [23JI increases only to about 3% 
for central and peripheral collisions in the calculated re- 
sults. For mid-central collisions there is a steady increase 
with the increase of transverse momentum to about 10% 
which is in line with the data. 

For the calculated pion elliptic flow at Ey^ = 40A GeV 
(Fig. I23p a quite similar behaviour as for protons is seen. 
At high p t the experimental results have large statistical 
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FIG. 20: Elliptic flow of protons in Pb+Pb collisions at i?i a b = 
20 A GeV and E la . h = 30 A GeV with p t < 2 GeV/c. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b — 5 — 9 fm for mid- 
central and b — 9 — 15 fm for peripheral collisions have been 
used. 



error bars. At low p t the experimental is not plagued by 
statistical and systematic errors and one observes good 
agreement between model and data. There is a increase 
of elliptic flow with increasing transverse momentum be- 
cause the larger the momentum of the particle the ear- 
lier it escapes the collision zone (53|. These high energy 
particles carry the signal of the very early stage of the 
collision where the coordinate space asymmetry is most 
pronounced (23| . 

For pion elliptic flow at 160 AGeV the picture is very 
similar to the plot of the proton flow at this energy. 

The predictions for the transverse momentum depen- 
dence of proton elliptic flow at £iab = 20A GeV and 
E\ a b = 30^4 GeV (Fig. ErJ)l have mostly the same shape 
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FIG. 21: Elliptic flow of pions in Pb+Pb collisions at £iab = 
2(L4 GeV and £ lab = 30A GeV with p t < 2 GeV/c. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b = 5 — 9 fm for mid- 
central and b = 9 — 15 fm for peripheral collisions have been 
used. 



as the calculations at 40^4 GeV. 

The elliptic flow of pions at 20A GeV and 30A GeV 
shown in Fig. [27] are also similar to the calculations for 
the pion elliptic flow at 40A GeV. But there is a difference 
for peripheral collisions. At the higher energy there is a 
lower value than for the lower energy. 



C. Excitation function 

The excitation function of charged particle elliptic flow 
is compared to data over a wide energy range (Fig. I28[) . 
i.e from E iah = 90 AMeV to = 200 GeV. The 

squeeze-out effect at low energies and the change to in- 
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FIG. 22: (Color online)Elliptic flow of protons in Pb+Pb col- 
lisions at S lab = 40 A GeV with \y\ < 1.1. UrQMD calcula- 
tions are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



plane emission at higher energies is nicely observed in 
the excitation function. The symbols indicate the data 
for charged particles from different experiments. Note 
however, that in the low energy regime there are only 
experimental data points for protons. For beam energies 
below 2 AGeV most of the charged particles are also pro- 
tons because there is not enough energy to produce many 
new particles. Going to higher energies the elliptic flow of 
pions and charged particles are very similar. The rapidity 
cut of \y\ < 0.1 has been used for the whole energy range 
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FIG. 23: (Color online) Elliptic flow of pions in Pb+Pb col- 
lisions at E^b = 40 A GeV with \y\ < 1.1. UrQMD calcula- 
tions are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 
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FIG. 24: (Color online)Elliptic flow of protons in Pb+Pb col- 
lisions at £i ab = 160A GeV with \y\ < 1.6. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of 6 < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



despite the fact that the data at higher energies is within 
\y\ < 0.5. This has been done to avoid too much changes 
in the parameters and this choice gives reasonable results 
over the whole energy range. We have checked that the 
results at higher energies are not affected by the choice 
of this narrower rapidity window. 

At low energies E\, eam ~ 0.1 — 6 A GeV the squeeze- 
out effect, i.e. the elliptic flow out-of-plane, is clearly 
seen in the data as well as in the calculations, especially 
when the mean field is considered. At such energies, 



it is well known that both the mean field and the two- 
body collision are equally important to reproduce quan- 
titatively the experimental results [HI, [6J, Hlj]. In this 
paper we adopt a hard equation of state with momen- 
tum dependence (HM-EoS) which was updated recently 
in UrQMD model [5l[. Meanwhile, the two-body scatter- 
ing in heavy ion collisions might be modified by the nu- 
clear medium. In order to consider (partly) the medium 
effect, the nucleon-nucleon elastic scattering cross sec- 
tions are modified to depend on the nuclear density, the 
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FIG. 25: (Color online) Elliptic flow of pions in Pb+Pb col- 
lisions at £i ab = 160/1 GeV with \y\ < 1.6. UrQMD calcu- 
lations are depicted with black lines. The symbols are NA49 
data from different analysis methods. The standard method 
(circles), cumulant method of order 2 (squares) and cumulant 
method of order 3 (triangles) are depicted. The 12.5% most 
central collisions are labeled as central, the centrality 12.5% 
-33.5% as mid-central and 33.5% -100% as peripheral. For 
the model calculations the corresponding impact parameters 
of b < 3.4 fm for central, b = 5 — 9 fm for mid-central and 
b = 9 — 15 fm for peripheral collisions have been used. 



isospin-asymmetry and the two-nucleon relative momen- 
tum, besides the center-of-mass energy of two nucleons. 
This treatment was investigated based on the relativistic 
Dirac-Brueckner-Hartree-Fock (DBHF) theory as well as 
the relativistic mean field (RMF) theory, please see Ref. 
[52| for details. Here we show the calculation results with 
the HM-EoS and with the DBHF-like medium modifica- 
tion on nuclcon-nulceon elastic cross sections (HMw) [13] . 

In the SPS regime the model calculations are quite 
in line with the data, especially with the NA49 results. 
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FIG. 26: Elliptic flow of protons in Pb+Pb collisions at 
E la . h = 20 A GeV and E lab = 30 A GeV with \y\ < 1. UrQMD 
calculations for 20 AGeV are depicted with solid lines while 
the results for 30 AGeV are depicted by dashed lines. Impact 
parameters of b < 3.4 fm for central, b — 5 — 9 fm for mid- 
central and b — 9 — 15 fm for peripheral collisions have been 
used. 



Above i^iab = 160 A GeV the calculation underestimates 
the elliptic flow. At the highest RHIC energy there are 
about 5% flow in the data while the model calculation 
provides only half of this value. This can be explained 
by assuming a lack of pressure in the transport model at 
these energies. 

In Fig. [29] the excitation function of elliptic flow of 
nucleons is presented. The rapidity cut of \y\ < 0.1 has 
been used because this is the appropriate one to compare 
with the data at lower energies. Due to this cut now also 
the calculation in the cascade mode (without nuclear po- 
tential) reaches negative values at low energies. With 
included potential the model is in line with experimen- 
tal data. The UrQMD result in the SPS energy range 
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FIG. 27: Elliptic flow of pions in Pb+Pb collisions at i?i a b = 
20A GeV and E la . h = 30A GeV with \y\ < 1. UrQMD cal- 
culations for 20 AGeV are depicted with solid lines while the 
results for 30 AGeV are depicted by dashed lines. Impact pa- 
rameters of b < 3.4 fm for central, b — 5 — 9 fm for mid-central 
and b = 9 — 15 fm for peripheral collisions have been used. 



lies in between the NA49 measurement of elliptic flow of 
protons. These values have been extracted from the dif- 
ferential plots of elliptic flow over rapidity for mid-central 
collisions discussed above in this paper. 

The observed proton flow V2 below ~ 5A GeV is 
smaller than zero, which corresponds to the squeeze-out 
predicted by hydrodynamics long ago 0, H, [H, [6?], IH, 
[691 ]. At higher energies, 10 — 160 A GeV, an increase of 
the flow vi to a maximum around £\ a b = 10A GeV fol- 
lowed by a decrease to about 2 % and a saturation is 
predicted from the string-hadronic transport model. In 
fact, the 158 A GeV data of the NA49 Collaboration sug- 
gest that a smooth increase proceeds between AGS and 
SPS. 
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FIG. 28: (Color online) The calculated energy excitation func- 
tion of elliptic flow of charged particles in Au+Au/Pb+Pb col- 
lisions in mid-central collisions (b=5-9 fm)with \y\ < 0.1 (full 
line). This curve is compared to data from different experi- 
ments for mid-central collisions. For E895 [H, HH and FOPI 
[5(| there is the elliptic flow of protons and for NA49 [Till i t 
is the elliptic flow of pions. For E877, CERES [13, IslEl, 
PHENIX [U, PHOBOS[Uj and STAR [f^there is data for 
the charged particle flow. The dotted line in the low energy 
regime depicts UrQMD calculations with the mean field [52] ]. 



The "collapse" of V2 (strong negative value) for pro- 
tons around midrapidity at AO A GeV is only pronounced 
in the standard method data. The UrQMD calculations, 
without a phase transition, show a robust 3% flow of 
protons. One cannot say anything about a clear underes- 
timation at high energies in this case because integrated 
proton flow data at RHIC is still not available. 



VII. SUMMARY 

We have compared UrQMD calculations to recent 
NA49 data. In general, a good agreement between data 
and calculation is found. There seem to be systematic 
uncertainties in the measurement method looking at the 
different results. For example for the directed flow data 
the effect of momentum conservation on the flow data 
can be seen. The slope around midrapidity of the rapid- 
ity distributions of proton directed flow becomes negative 
around E\ a \> — AO A GeV. This cannot be reproduced by 
the transport model calculations. The excitation func- 
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tion of elliptic flow shows strong negative flow at low 
energies - the "squeeze-out" -effect - which can quantita- 
tively only be reproduced by including a nuclear potential 
in the calculation. At high energies we observed an un- 
derestimation of the elliptic flow of charged particles in 
the present model. This can possibly be explained by 
assuming a lack of pressure in the early stage of the col- 
lisions at high energies. It will be very interesting to see 
what happens in the lower energy regime (20/30 A GeV) 
when high quality CBM-FAIR data becomes accessible. 
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